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ABSTRACT 

We present an analysis of the highly variable accreting X-ray pulsar 3A 1954+319 using 2005-2009 moni- 
toring data obtained with INTEGRAL and Swift. This considerably extends the pulse period history and covers 
flaring episodes in 2005 and 2008. In 2006 the source was identified as one of only a few known symbiotic 
X-ray binaries (SyXBs), i.e., systems composed of a neutron star accreting from the inhomogeneous medium 
around an M-giant star. The extremely long pulse period of ~5.3 hr is directly visible in the 2008 INTEGRAL- 
ISGRI outburst light curve. The pulse profile is double peaked and generally not significantly energy dependent 
although there is an indication of possible softening during the main pulse. During the outburst a strong spin- 
up of —1.8 x 10~ 4 hrhr 1 occurred. Between 2005 and 2008 a long-term spin-down trend of 2.1 x 10~ 5 hrhr 1 
was observed for the first time for this source. The 3-80 keV pulse peak spectrum of 3 A 1954+319 during the 
2008 flare could be well described by a thermal Comptonization model. We interpret the results within the 
framework of a recently developed quasi-spherical accretion model for SyXBs. 

Subject headings: binaries: symbiotic — stars: individual (3A 1954+319) — stars: neutron — X-rays: binaries 



1. INTRODUCTION 

The X-ray source 3 A 1954+3 19 was detected in the Cygnus 
region in surveys by Uhuru, Ariel V, EXOSAT, and ROSAT 
(Forman et al. 1978; Warwick et al. 1981, 1988; Voges et al. 
1999). Pointed observations with EXOSAT (Cook et al. 1985) 
and Ginga (Tweedy et al. 1989) showed a hard X-ray spec- 
trum as well as intensity variations by an order of magnitude 
on timescales of minutes. This led to the suggestion that the 
system might be a High Mass X-ray Binary (HMXB). Only 
when Masetti et al. (2006) identified the companion as an 
M4-M5 III star at a distance of <1.7kpc and Corbet et al. 
(2006) discovered a ~5 hr pulse period in early Swift-BAT 
data, did it become clear that 3A 1954+319 is a Symbiotic X- 
ray Binary (SyXB). SyXBs constitute a small group 7 of per- 
sistent Low Mass X-ray Binaries (LMXBs) in which a neu- 
tron star is orbiting in the inhomogeneous medium around an 
M-type giant star 8 . 

SyXBs typically have wide orbits, e.g., the prototype 
GX1+4 has an orbital period of ~ 1161 days (Hinkle et al. 
2006). Their X-ray emission is therefore due to wind accre- 
tion, a process not well investigated for late type donors. Two 
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broadband spectral studies support the SyXB interpretation 
for 3A 1954+319: Mattana et al. (2006) modeled the non- 
simultaneous BeppoSAX and INTEGRAL spectrum with a 
highly absorbed (Nu ~ 10 23 cirr 2 ) cutoff power law with a 
photon index of 1.1 and a folding energy of 15keV and a 
weak Fe Ka line. These authors also confirmed the detec- 
tion of the ~5 hr period in BAT and ISGRI data. In a study of 
archival data spanning absorbed 2-10 keV luminosities from 
3.4xl0 34 ergs-' to 1.8xl0 35 ergs^ Masetti et al. (2007) con- 
firmed the empirical spectral description for the >2 ke V spec- 
trum and determined a best fit using thermal Comptonization, 
modified by complex absorption (ionized plus partially cov- 
ering neutral absorption). They also modeled a <2keV soft 
excess with a ~50 eV hot plasma and interpreted it, together 
with the two-zone absorption, as being due to a diffuse, partly 
ionized cloud of material around the neutron star. 

The ~5 hr period is the only period known for this system. 
As Corbet et al. (2006, 2008) argued (i) the period value itself 
is inconsistent with being the orbital period of an M-giant, 
(ii) the large period decline of (2.6 ± 0.2) x 10 _5 hrhr _1 ob- 
served over the first year of BAT data cannot be due to or- 
bital Doppler modulation, and (iii) the period change is also 
too large to be supported by a white dwarf accretor. Inter- 
preting the ~5 hr period as a neutron star spin period makes 
3 A 1954+319 the slowest rotating neutron star in an X-ray bi- 
nary currently known and one of the slowest pulsars in general 
- with only the 6.67 hr pulsar in the supernova remnant RCW 
103 showing a larger value (De Luca et al. 2006). Corbet 
et al. (2008) also noted that if the neutron star rotated close to 
its equilbrium period for disk accretion, as usually assumed in 
accreting sources, the long period would imply a neutron star 
magnetic field of ~10 15 G rather than of ~10 12 G. The latter 
value is more commonly observed for accreting pulsars and 
is consistent with a spin-down origin of the long period (Mat- 
tana et al. 2006). From the Swift-BAT and RXTE-ASM data 
available at the time it was not clear whether the strong spin- 
up observed was associated with the long flaring episode in 
2005 and/or whether the source might show a spin-down at 
lower fluxes, as it is the case for GX 1+4 (Corbet et al. 2008). 
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Figure 1. From top to bottom: INTEGRAL observations of 3A 1954+319 
with an offset angle <10° and INTEGRAL -ISGRI detections (blue and brown 
tickmarks, respectively). The long-term light curve shown was obtained by 
Swift-BAT in the 15-50 keV range and has been rebinned to a resolution 
of 5 days. The lower part of the figure shows the pulse period evolution as 
determined by BAT (black) and ISGRI (red). 

(A color version of this figure is available in the online journal.) 

In this Letter we considerably extend the pulse period his- 
tory for 3A 1954+319 and present a timing and spectral anal- 
ysis of a flaring episode in 2008 (Fiirst et al. 2011 presented 
a preliminary analysis of these data). Section 2 describes the 
observations and data reduction. Section 3 reports the results, 
including long-term and high time resolution light curves, en- 
ergy resolved pulse profiles, the pulse period evolution, and 
the broadband spectrum. In Section 4 the results are summa- 
rized and discussed. 

2. OBSERVATIONS AND DATA REDUCTION 

The upper part of Figure 1 shows the 2005-2009 15- 
50keV light curve of 3A 1954+3 19 9 observed with Swift- 
BAT (Barthelmy et al. 2005). Months long flaring episodes 
are apparent, especially in 2005 and in 2008. The former 
includes most of the time range analyzed by Corbet et al. 
(2008), -MJD 53330-53680. The short tickmarks above the 
BAT light curves indicate 1163 INTEGRAL (Winkler et al. 
2003) pointings ("science windows", ~2 ks exposures) during 
which 3 A 1954+319 was within the field of view of the ISGRI 
(Lebrun et al. 2003) of the IBIS instrument, with a pointing 
offset < 10°. 

Since its launch in 2002 INTEGRAL has performed sev- 
eral extensive monitoring campaigns of the Cygnus region 
(Pottschmidt et al. 2003; Cadolle Bel et al. 2006; Martin et al. 
2009; Laurent et al. 201 1; Williams et al. 201 1). A Key Pro- 
gram (KP) centered on the black hole X-ray binary Cygnus X- 
1, located 3°. 15 from 3A 1954+319, has been in place since 
2008 with annual exposures of several lOOks (Grinberg et al. 
2011). Coincidentally these KP observations covered about 
three weeks of the 2008 flare of 3A 1954+319 in unprece- 
dented detail (Figure 1). 

The 201 science windows for which 3A 1954+319 was de- 
tected (OSA DETSIG >6 in the 20-100 keV science window 
images) are indicated by long tickmarks in Figure 1 . For these 
pointings the ii_light tool in version 7 of the Offline Scien- 
tific Analysis (OSA; Courvoisier et al. 2003) 10 was applied to 

9 Daily light curve from http://heasarc.gsfc.nasa.gov/docs/ 
swift/results/transients/4U1954p31/. 

10 The ii.light tool was not present in OSA 8 and at the time of writing 
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Figure 2. Upper panel: INTEGRAL -ISGRI 20-100 keV light curve of 
3A 1954+319 observed during the flaring episode in 2008. Each data point 
represents one science window (~2ks). The insets show two close-ups with 
the full 100 s resolution light curve (gray), the count rates per science window 
(magenta), and the average pulse profile (red histogram). Lower panel: The 
three pulse period measurements obtained with ISGRI during the flare. 

(A color version of this figure is available in the online journal.) 

produce ISGRI light curves of 3 A 1954+319 with a time res- 
olution of 100 s in the energy bands 20-40 keV, 40-100 keV, 
and 20-100 keV. Together with the BAT long-term light curve 
the 20-100 keV ISGRI light curves were used to determine 
the pulse period evolution (Section 3.3). 

In addition, a more detailed analysis was performed for the 
flare in 2008. This dataset included science windows from 
satellite revolutions 739, 741-746, 756, and 758 (one revolu- 
tion takes about three days). Pulse profiles in the three energy 
bands were created (Section 3.2). In order to maximize the 
signal to noise ratio (S/N) all science windows in the phase 
range of 0.45-0.85, i.e., associated with the main pulse peak, 
were selected. OSA 9 was used to extract average ISGRI 
and Joint European X-ray Monitor (JEM-X; Lund et al. 2003) 
spectra, for source offset angles <10° and <3°, respectively. 
The ISGRI spectrum had an exposure of 85 ks, was created 
by averaging spectra from individual science windows, and 
was modeled in the 20-80 keV range. The JEM-X spectrum 
had an exposure of 9.8 ks, was extracted from mosaic images 
using mosaic_spec and was modeled in the 3-30keV range 
(Section 3.4). Response and auxiliary response files were se- 
lected or created following the analysis documentation 11 . 

3. RESULTS 
3.1. Light Curves 

The long-term BAT light curve shown in Figure 1 demon- 
strates the irregular flaring of 3A 1954+319 on timescales 
of months, with the 2008 November flare having been one 
of the brightest since the start of the BAT monitoring. The 
INTEGRAL observations cover the second half of the flar- 
ing episode. The resulting 20-100 keV ISGRI light curve is 
shown in the upper panel of Figure 2. The 20-100 keV flux 
during the outburst varied by a factor of ~20 with an average 
of ~40mCrab and a peak value of ~130mCrab. 

The insets of Figure 2 show close-ups of two randomly se- 
lected parts of the ISGRI outburst light curve with a resolu- 
tion of 100 s. Individual pulses are directly observed in the 
light curve for the first time. Also shown are repetitions of the 

was still undergoing evaluation for OSA 9. 

1 1 See http://www.isdc. unige.ch/integral/analysis#Documentation 
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Figure 3. Upper panel: ISGRI pulse profiles for the flaring episode in 2008, 
in the energy ranges of 20-40 keV (black) and 40-100 keV (red, multiplied 
by 3 for better visibility). Middle panel: Hardness ratio obtained by divid- 
ing the 40-100 keV by the 20-40 keV profile. The dashed line indicates the 
mean hardness. Lower panel: Deviation of the hardness ratio from the mean 
hardness in units of cr. The dashed line indicates no deviation. 
(A color version of this figure is available in the online journal.) 

P-corrected average pulse profile obtained by folding the high 
resolution light curve on the pulse ephemeris determined from 
the outburst data (Section 3.3). Comparing the profiles with 
the high time resolution light curve and the average science 
window count rates demonstrates general consistency but al- 
lows for moderate pulse-to-pulse variations which are com- 
mon in accreting X-ray pulsars (Klochkov et al. 201 1). 

3.2. Pulse Profiles 

Pulse profiles in the energy ranges of 2CMM) keV and 40- 
lOOkeV were obtained by folding the energy resolved high 
time resolution ISGRI outburst light curves on the pulse 
ephemeris determined from the 2008 outburst data (Sec- 
tion 3.3), see upper panel of Figure 3. These are the highest 
quality pulse profiles available for the source to date. They 
clearly show a double peaked structure, in contrast to the sin- 
gle peaked pulse profile displayed by the prototype SyXB 
GX 1+4 (Ferrigno et al. 2007). While the possible presence 
of the secondary peak for 3 A 1954+319 was indicated in the 
<50keV BAT profiles presented by Corbet et al. (2008), it 
could not be detected in the 50-100 keV band, possibly due 
to the comparatively smaller S/N. The lower two panels of 
Figure 3 show the hardness ratio between the profiles in the 
two energy bands and its deviation from the average ratio. No 
significant energy dependence was detected with exception of 
a possible moderate softening during the brightest part of the 
main pulse in a narrow phase range (~ 0.60 - 0.65). 

3.3. Pulse Period Evolution 

Using the epoch folding technique (Schwarzenberg-Czerny 
1989) pulse period values were determined for three equally 
long parts of the 2008 outburst observations with ISGRI. The 
results are shown in the lower panel of Figure 2. Uncertain- 
ties were calculated according to the Monte Carlo method de- 
scribed by Davies (1990): for every segment 10 4 light curves 
with the same sampling and variance as seen in the observa- 
tional data were simulated, based on the segment's average 
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Figure 4. Upper panel: JEM-X (3-30 keV, brown) and ISGRI (20-80 keV, 
blue) counts spectra for the flaring episode in 2008. In order to obtain a better 
S/N only data from the main pulse, i.e., phases 0.45-0.85, were included. The 
best fit description consisted of an absorbed thermal Comptonization model 
and is also displayed. Lower panel: Best fit residuals. 

(A color version of this figure is available in the online journal.) 

pulse profile modified by Gaussian noise. A successful pe- 
riod search was performed for every simulated light curve and 
the widths of the emerging distributions were used as uncer- 
tainties of the pulse period measurements. The three ISGRI 
pulse period values for the flare evolved from 5.336 ± 0.003 hr 
over 5.308 + 0.003 hr to 5.264 ± 0.006 hr, i.e., a strong spin- 
up became apparent. The pulse ephemeris over the 2008 
flare was determinded to T = MJD54782.6897, P(T ) = 
5.3060 ± 0.0007 hr, and P(T ) = (-1.81 ±0.17)x lO^hrhr" 1 
(uncertainties are given on a Icr confidence level). 

In order to put the ISGRI period values into perspective 
the pulse period history from the long-term BAT light curve 
was updated by performing local period determinations for 
20 days long segments. The lower part of Figure 1 contains all 
successful BAT period measurements, the three ISGRI flare 
measurements, as well as three additional ISGRI measure- 
ments that could typically be obtained during times of denser 
sampling and/or elevated count rates 12 . The long-term results 
are the following: (i) the spin-up phase during the flaring ac- 
tivity in 2005, analyzed by Corbet et al. (2008) and Mattana 
et al. (2006), was reproduced, (ii) it was followed by a long 
spin-down trend between 2005 and 2008, characterized by 
P ~ 2.1 x 10- 5 hrhr _1 , (iii) the BAT and ISGRI pulse pe- 
riod values are in excellent agreement, especially during the 
strong spin-up in 2008, (iv) the spin-down trend resumed in 
2009, possibly slowed down by continued moderate flaring. 
Note that while the possibility of a beginning spin-down at 
the end of 2005 was mentioned by Corbet et al. (2008), a 
spin-down has now been clearly observed for the first time 
for 3A 1954+319. Also note that the spin-up in 2008 was an 
order of magnitude larger than the one in 2005. 

3.4. Broadband Spectrum 

As described above the 3-80 keV spectrum of 
3A 1954+319 was determined for the phase range asso- 
ciated with the peak of the main pulse. This approach was 
chosen because the source is close to the detection limit 
during dimmer pulse phases. The spectrum could be well 

12 Following a conservative approach only peaks which were visually 
clearly apparent in the epoch folding statistics were selected. 
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described by an absorbed thermal Comptonizaton model (Fig- 
ure 4), resulting in x 2 eA = 1.3 for 32 degrees of freedom. The 
best fit parameters obtained were A^h = 3.9^ 2 x 10 22 cm~ 2 , 
kT = 0.9!^keV, t — 4.2 + 0.7, and kT e = 7.5^ keV. 
A flux cross calibration factor of 0.84+9J 2 was applied to 
model the JEM-X data. The residuals mainly reflect known 
ISGRI calibration uncertainties (Grinberg et al. 2011). 
All fit parameter uncertainties, including those of the flux 
measurements below, are given at the 90% level for one 
interesting parameter. A broken power law fit provided a 
slightly better description but the improvement is due to 
masking cross-calibration issues (x 2 ~ 1.2). A cutoff power 
law model as applied by Mattana et al. (2006) resulted in a 
worse fit (y 2 ~ 1.9). 

The absorbed 2-10keV and 10-100 keV fluxes were 7.2+° * 
and 12.1*Q 3 ! XlO _10 ergcrnr 2 s _1 , respectively (ISGRI normal- 
ization). The unabsorbed fluxes were 9.1*|q and 12.2^3 x 
10~ 10 erg cirT 2 s _1 , corresponding to (pulse peak) luminosities 
of 3.1 and 4.2 x 10 35 erg s _1 for a distance of 1.7 kpc. Since the 
assumed distance is an upper limit (Masetti et al. 2006), these 
luminosities are upper limits as well. They are comparable to 
the highest luminosities reported by Masetti et al. (2007). 

4. DISCUSSION 

A possible scenario for the development of pulse periods 
longer than a few hundred seconds has been proposed for 
persistently bright HMXBs by Ikhsanov (2007), where accre- 
tion proceeds spherically and a strong spin-down happened 
in a previous accretion epoch (subsonic propeller regime). 
The author argues that the 2.7 hr pulse period of the HMXB 
2S 01 14+650 (Farrell et al. 2008) can thus be explained with- 
out the need for an unusually high (10 15 G) magnetic field. 

For the LMXB 3 A 1954+319 the long pulse period is more 
easily reconcilable with the system's life time (Mattana et al. 
2006). Shakura et al. (2011) recently developed an accre- 
tion model for SyXBs that can not only account for long spin 
periods but also provides a mechanism allowing for quasi- 
spherical accretion. In this model a subsonic settling regime 
occurs for X-ray luminosities below ~3 x 10 36 erg s -1 , i.e., for 
luminosities consistent with those observed for 3 A 1954+319. 
A shell of hot material forms around the magnetosphere 
which mediates the transfer of angular momentum to/from the 
neutron star by advection and viscous stress. The accretion 
rate is determined by the ability of the plasma to enter the 
magnetosphere. The equilibrium pulse period in this case is: 

Peg / B \ n/U ( P oMt \( L x \ 4/ "/ y wlnd V* 
10 3 s \10 12 G/ I lOd/UO^ergs- 1 / \10 3 km S - 1 / 

(1) 

where P 01 -bit is the binary period and v w i nl j is the stellar wind 
velocity with respect to the neutron star. Reporting on popu- 
lation synthesis simulations for SyXBs based on this model, 
Postnov et al. (2011) showed that the 5.3 hr pulse period of 
3 A 1954+319 can be well reproduced. 

3A 1954+319 seems to show long spin-down episodes be- 
tween major flares (P ~ 2.1 x lO^hrhr -1 ) which are remi- 
niscent of the spin-down displayed by GX 1 +4 since the early 
1980s (P ~ lO^hrhr 1 ; Gonzalez-Galan et al. 2011). Sev- 
eral models beyond equilibrium disk accretion have been pro- 
posed to explain the spin-down of GX 1 +4 without invoking 
unusually high magnetic fields, for example the presence of a 
counterrotating disk (Nelson et al. 1997) or accretion of fall- 



back material expelled during the propeller phase (Perna et al. 
2006). Only the quasi-spherical accretion model, however, re- 
produces the correct sign and magnitude of the negative corre- 
lation between spin frequency change v and X-ray flux Fx ob- 
served during spin-down in this source (~v <x F^! 7 ; Gonzalez- 
Galan et al. 201 1). For 3A 1954+319 such a detailed study of 
the v-Fx-relationship is difficult, especially during spin-down 
since the low flux allowed for only a few v measurements 
(Figure 1). It is beyond the scope of this Letter. 

On longer time scales we observe torque reversals and a 
positive v-Fx-correlation between low (spin-down) and high 
flux (spin-up) episodes (Figure 1). This is consistent with the 
behavior predicted by the quasi-spherical accretion model for 
higher accretion rates within the settling regime (see Figure 1 
of Shakura et al. 2011). 

The strong spin-up in 2008 translates to P/P = -0.9 x 
10~ 8 s~'. While still high, the absolute value is of the 
same order of magnitude as the spin-down related P/P of 
3.1 x 10~ 8 s- 1 observed for the SyXB IGR J16358-4724 (Pa- 
tel et al. 2007). The X-ray luminosity required to sustain such 
a spin-up in 3A 1954+319 in the equilibrium disk accretion 
case would be ~5 x 10 36 ergs~' (Joss & Rappaport 1984), 
whereas we obtained an upper limit for the pulse peak flux 
value of ~7.4 x 10 35 erg s _1 , again arguing against disk accre- 
tion. 

This conclusion does not change when the possible contri- 
bution of a spin-up due to orbital motion is considered: Ac- 
cording to Dumm et al. (1998) the measured maximum mass 
for M giants is ~ 3.5 M and the median stellar radii deter- 
mined for M4 III and M5 III stars are 103 R Q and 12Ofi , 
respectively. For neutron star orbits outside of the M giant, 
the orbital period then has to be > 80 days and > 100 days, 
translating into a maximum relative change of the pulse pe- 
riod over the orbit of < 4 x 10~ 4 in both cases. The lowest 
relative uncertainty range of our pulse period measurements 
is ~ 10~ 3 , considerably bigger than any realistic orbital effect, 
which will be even smaller than the value above due to a typi- 
cally lower stellar mass and wider orbit. Therefore the orbital 
influence is negligible for the 2008 spin-up measurement. 

The best fit to the broadband spectrum of the 2008 flare de- 
scribes an optically thick Compton plasma (r = 4.2 ± 0.7, 
lcT e = 7.5+Q4keV) with parameters qualitatively consistent 
with the results of Masetti et al. (2007). As these authors 
state the parameters are similar to those commonly seen in 
LMXBs with an accreting neutron star companion. Most of 
those sources do not show pulsations, however, and their ac- 
cretion geometry is most likely different from that of pul- 
sars. Other SyXBs show similar spectra, especially GX 1+4 
(t = 6.80 ± 0.15, kT e = 13.1 ± 0.2 keV; Ferrigno et al. 2007). 
Furthermore, optically thick Comptonization has also been 
used to describe the emission from the accretion columns of 
accreting pulsars in HMXBs, e.g., for the cyclotron line soure 
1A 1118-61 (Suchy et al. 2011). We tentatively propose a 
similar origin close to the neutron star surface for the broad- 
band X-ray emission of 3A 1954+319. 
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